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Introduction 

Electric  arc  heaters  produce  high-temperature  gas  which  can  be  used  for  a 
number  of  purposes,  including  the  generation  of  thrust  for  spacecraft  propul¬ 
sion,  the  synthesis  of  diamond  and  boron  nitride  films  for  industrial  tool 
hardening  and  electrical  insulation,  respectively,  and  the  generation  of  high- 
enthalpy,  hypersonic  air  streams  for  the  investigation  of  processes  associated 
with  ultra-high  velocity  atmospheric  flight.  As  ubiquitous  as  electrical  arcs 
are,  understanding  of  them  is  limited.  As  such,  attempts  at  extending  the  op¬ 
erating  envelope  of  electrical  arc  devices  (e.g.  by  increasing  pressure)  through 
empiricism  alone  can  lead  to  unexpected  and  sometimes  disastrous  results. 

The  work  described  here  is  a  combined  numerical/experimental  study  of 
arcjets.  The  focus  of  the  experimenal  work  has  been  a  parametric  study 
of  the  onset  of  instability  in  arc  thruster  devices.  The  focus  of  the  numer¬ 
ical  work  has  been  the  development  of  numerical  methods  for  solving  the 
magnetohydrodynamics  (MHD)  equations  in  complex  geometries. 
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Summary  of  Experimental  Work 

High-enthalpy  arc-tunnels  enter  an  unstable  operating  regime  when  the  work¬ 
ing  fluid  mass  flow  rate  is  maintained  at  an  incompatible  value.  A  similar 
phenomenon  is  manifested  in  arc-producing  propulsion  systems  (e.g.,  arcjets) 
in  a  process  known  as  starvation,  whereby  the  propellant  flow  rate  is  too  low 
to  adequately  satisfy  particle  flux  requirements  of  a  stable  arc.  In  order  to 
determine  how  this  phenomenon  and  others  seen  in  arcjets  may  be  used  to  un¬ 
derstand  arc-tunnel  instabilities,  the  discharge  current  and  voltage  of  a  1  kW 
arcjet  were  recorded  as  a  function  of  mziss  flow  rate  by  a  digital  oscilloscope. 
This  Wcis  done  to  determine  how  propellant  flow  rate  affects  the  stability  of 
the  arcjet  dischajge  and  to  see  how  this  information  may  be  used  in  the  de¬ 
velopment  of  a  computer  model.  These  measurements  were  performed  with 
hydrogen  as  the  propellant.  In  addition,  near-field  electron  temperature  and 
number  density  measurements  were  made  via  triple-Langmuir  probe  for  com¬ 
parison  against  code  predictions.  A  1  kW  arcjet  (Fig.  1)  that  was  supplied 
by  the  NASA  Lewis  Research  Center  (LeRC)  was  used  for  this  study.  The 
arcjet  was  first  developed  as  a  by-product  of  arc  [high-enthalpy]  tunels  used 
in  the  early  space  age  to  develop  heat  shields.  An  arcjet  uses  an  arc  to  heat 
propellant  (e.g.,  hydrogen)  which  is  subsequently  expanded  through  a  nozzle. 
The  arcjet  used  in  this  study  features  a  2%-thoriated  tungsten  cathode  and 
a  nozzle  (also  of  2%-thoriated  tungsten)  that  serves  as  the  anode.  The  arcjet 
has  a  0.51-mm-diameter  by  0.25-mm-long  constrictor,  a  30-degree  half-angle 
converging  nozzle  section  upstream  of  the  constrictor,  and  a  20-degree  half¬ 
angle  diverging  section.  The  exit  diameter  of  the  nozzle  is  9.5  mm,  giving 
the  expansion  section  an  area  ratio  of  350.  The  electrode  gap  spacing  is  0.51 
mm  and  the  outer  housing  of  the  device  is  constructed  of  titaniated  zirco- 
niated  molybdenum.  The  nominal  exhaust  velocity  is  approximately  8,000 
m/s  with  hydrogen.  All  experiments  reported  were  performed  in  the  Plasma- 
dynamics  and  Electric  Propulsion  Laboratory  6  m  by  9  m  vacuum  chamber 
(Fig.  2).  At  the  time  of  these  tests,  the  chamber  was  pumped  with  six 
81-cm-diamter  oil  diffusion  pumps,  two  blowers  and  four  mechanical  pumps, 
providing  a  pumping  speed  of  300,000  1/s  on  hydrogen.  Background  chamber 
pressure  was  maintained  to  less  than  1.5  x  10-4  Torr  during  arcjet  operation. 
Since  then,  four  CVI  TM1200  reentrant  cryopumps  has  been  installed  on 
the  facility.  Arcjet  power  was  provided  by  a  1800  VV  Sorenson  power  supply 
that  was  conditioned  by  a  NASA  LeRC  power  processing  unit  (PPU).  Arcjet 
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voltage  was  measured  with  Tektronix  P6007  100:1  voltage  probes  clamped 
to  the  electrode  leads  on  the  vacuum  chamber  bulkhead.  The  voltage  probe 
signals  were  collected  by  a  Tektronix  AM501  operational  amplifier.  Arcjet 
current  is  monitored  with  a  Tektronix  A6303  current  sensor  powered  by  a 
Tektronix  AM503  current  probe  amplifier.  Arcjet  voltage  and  current  data 
were  monitored  and  recorded  by  a  Tektronix  TDS  540  digitizing  oscilloscope, 
and  stored  on  a  computer  via  a  GBIP  interface  and  Lab  VIEW  data  acqui¬ 
sition  system.  The  arcjet  was  operated  at  hydrogen  flow  rates  of  up  to  10 
SLM  (15  mg/s).  Propellant  was  supplied  to  the  arcjet  from  compressed  gas 
bottles  through  stainless-steel  feed  lines.  Propellant  flow  was  controlled  and 
monitored  with  an  MKS  1159B  mass  flow  controller  specially  calibrated  for 
light  gases.  A  triple  probe  (Fig.  3)  was  used  to  measure  ne  and  Te  near  the 
exit  plane  of  the  arcjet  while  running  at  10  A,  79  V,  and  with  a  flow  of  10 
SLM.  The  probe  was  comprised  of  three  individual  tungsten  wire  electrodes, 
parallel  to  each  other,  0.23  mm  in  diameter  with  an  exposed  length  of  5.0 
mm.  Each  wire  was  mounted  in  round  single-bore  alumina  tubes  which  were 
held  together  with  a  high-temperature  ceramic  adhesive.  The  spacing  be¬ 
tween  the  electrodes  was  1  mm.  Two  Duracell  alkaline  batteries  were  used 
to  supply  Vd3.  Vd2,  Vd3.  These  voltages  as  well  as  the  current  shunt  volt¬ 
age,  Vsh,  were  directly  fed  to  the  Lab  VIEW  data  acquisition  system.  The 
probe  was  mounted  on  a  large  positioning  system  which  contains  two  rotary 
platforms  on  a  6  ft-long  radial  stage  that  is  mounted  on  a  3  ft  travel  axial 
stage.  The  probe  was  not  rotated  to  be  aligned  with  the  flow  due  to  time 
limitation.  The  errors  due  to  probe  misalignment  with  the  flow  will  affect  the 
radial  measurements.  This  misalignment  may  result  in  artifically  high  Te  and 
ne  in  some  cases.  The  positioning  system  has  an  absolute  position  accuracy 
of  0.15  mm.  Experimental  Results  Discharge  Characteristics  The  discharge 
current  and  discharge  voltage  were  recorded  for  each  propellant  flow  rate  us¬ 
ing  a  digital  oscilloscope.  The  data  are  shown  in  Figures  4  through  7.  The 
horizontal  scale  (time  scale)  of  the  oscilloscope  traces  was  chosen  to  be  20 
s  per  division.  This  scale  was  selected  to  show  clearly  both  the  structures 
of  the  signals  and  the  changes  in  those  structures  with  varying  flow  rates 
at  the  same  time.  The  oscillation  of  about  11  kHz  in  the  discharge  current 
and  discharge  voltage  is  caused  by  the  oscillation  in  the  output  of  the  PPU. 
The  discharge  current  remains  constant  at  10  A  while  the  discharge  voltage 
decreases  as  the  flow  rate  decreases.  This  is  due  to  the  fact  that  the  PPU  is  a 
constant-current  power  supply.  The  oscilloscope  traces  clearly  indicate  that 
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the  structures  in  the  discharge  current  and  discharge  voltage  signals  become 
"rattled”  as  the  flow  rate  decreases  and  the  arc  becomes  starved.  It  is  not 
shown  here  that  as  the  flow  rate  decreased,  the  repeatability  of  those  signals 
suffered  as  well.  The  discharge  current  signals  were  less  affected  by  the  de¬ 
crease  of  the  flow  rate  than  the  discharge  voltage  signals.  This  is  most  likely 
due  to  the  constant-current  mode  of  the  PPU.  Although  there  is  increased 
hash  at  the  lower  flow  rates,  the  peak-to-peak  values  of  the  discharge  current 
and  discharge  voltage  signals  do  not  vary  significantly  while  the  flow  rate  is 
changed.  This  may  be  due  to  the  PPU  as  well.  The  oscilloscope  traces  also 
shows  that  there  is  a  certain  flow  rateor  a  range  around  a  certain  flow  rateat 
which  the  discharge  switches  to  a  "rattled  mode.”  In  this  study,  it  was  found 
to  be  at  about  4.75  SLM  of  H2  flow.  This  result  may  be  useful  when  one  tries 
to  save  as  much  propellant  as  possible  without  destabilizing  the  discharge. 
This  result  may  also  be  useful  in  the  development  of  a  computer  code  for 
arcjet  operation.  Plume  Measurements  Axial  and  radial  profiles  of  the  arcjet 
plume  while  the  arcjet  was  running  at  10  A/79  V  with  a  H2  flow  rate  of  10 
SLM  were  generated.  Axial  profile  measurements  of  ne  made  with  a  triple 
probe  along  the  axis  of  the  arcjet  from  5  mm  to  300  mm  downstream  of  the 
exit  plane  were  presented  in  Fig.  8.  These  measurements  show  the  electron 
temperature  and  number  density  to  vary  between  0.25  and  0.4  eV  (2900-4600 
K)  and  1011  to  2  x  1012  cm-3,  respectively.  Radial  profile  measurements  of 
ne  and  Te  made  with  the  same  triple  probe  at  5  mm  downstream  of  the  exit 
plane  were  presented  in  Fig.  9.  Figure  8  shows  the  comparison  of  centerline 
ne  results  from  this  study  and  a  previous  study.  1  The  two  studies  match  well 
beyond  40  mm  from  the  exit  plane.  Electron  number  densities  measured  in 
this  study,  however,  is  about  half  of  those  of  the  previous  work  for  closer  dis¬ 
tances  from  the  exit  plane.  The  power  level  of  the  arcjet  during  the  present 
experiments  was  0.79  kW,  which  is  about  20%  less  than  the  power  level  of 
1  kW  of  Reference  1.  Thus,  the  particle  flux  was  less  in  this  study  than  in 
the  other  study.  This  difference  in  particle  flux  due  to  lower  engine  power  is 
small  far  from  the  exit  plane  because  the  flux  of  particles  there  is  mostly  due 
to  diffusive  processes.  Fig.  9  shows  the  radial  profile  of  ne  and  Te.  It  seems 
that  there  are  two  distinct  regions  separated  at  the  radial  position  of  8  mm 
(more  apparent  in  the  Te  profile).  When  the  plume  was  observed  visually, 
there  were  two  distinct  regions  in  the  plume  in  a  ”co-annular  cone.”  The 
inner  cone  was  red,  and  the  surrounding  region  was  light  blue.  It  should  be 
noted  that  this  two-color  plume  was  not  observed  when  the  thruster  ran  at 
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1  kW  of  power.  However,  there  might  still  have  been  the  similar  structure 
in  the  plume,  though  not  readily  visible.  Conclusions  The  discharge  cur¬ 
rent/voltage  measurements  at  different  flow  rates  suggest  that  the  discharge 
becomes  unstable  as  the  propellant  flow  rate  decreases,  and  that  there  exists 
a  certain  flow  rate,  or  a  small  range  around  a  certain  flow  rate,  at  which 
the  discharge  switches  into  a  destabilized  mode.  The  plume  study  suggests 
that  there  is  a  co-annular-cone-shaped  structure  in  the  arcjet  plume.  The 
study  also  indicates  that  20  %  reduction  in  the  engine  power  results  in  50  % 
decrease  in  the  particle  flux  near  the  exit  plane. 
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Summary  of  Numerical  Work 


The  numerical  modelling  focused  on  solution  of  the  ideal  magnetohydrody¬ 
namic  equations.  The  equations  to  be  solved  are  those  for  conservation  of 
mass,  momentum,  magnetic  field  and  energy.  They  are  solved  in  the  form' 
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which,  with  the  exception  of  the  source  term  on  the  right-hand  side  of  the 
equation,  is  the  divergence  form  typically  in  solving  systems  of  conservation 
laws.  The  source  term  is  present  so  that  any  local  deviations  from  V  •  B  are 
handled  correctly.  In  this  form,  the  quantity  (V  •B)/p  is  treated  as  a  passive 
scalar;  any  (V  •  B)//!)  that  is  created  is  passively  convected.  This  form  of 
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the  equations  is  the  symmetrizable  form,  and  has  the  advantage  of  Gallilean 
invariance  and  the  ability  to  derive  an  entropy  conservation  law  by  a  linear 
combination  of  the  remaining  equations.  Details  of  the  governing  equations 
and  the  solution  scheme  are  given  in  the  attached  paper,  which  has  been 
submitted  to  Journal  of  Computational  Physics. 

Equation  1  is  solved  by  a  Roe  scheme  (Reference  1),  an  HLLE  scheme 
(Referenced)  or  a  Boltzmann  scheme  (Reference  5)  derived  from  the  eigensys- 
tem  of  the  quasilinear  form  of  Equation  1.  Time-stepping  is  done  explicitly, 
and  the  entire  algorithm  has  been  tied  into  a  solution-adaptive  code.  Rep¬ 
resentative  results  are  shown  in  Figures  10-13.  Figure  10  shows  the  grid 
used  in  the  set  of  calculations;  the  geometry  is  a  simple  arc  tunnel.  The 
grid  is  adaptively  refined  automatically  so  as  to  resolve  the  anode,  cathode, 
and  the  arc.  Figures  11-13  show  the  current  density,  acoustic  Mach  number 
and  plasma  density  in  the  arc  tunnel  for  the  three  solvers.  The  Boltzmann 
scheme  is  the  most  diffusive,  but  also  the  least  expensive  computationally. 
The  HLLL  scheme  is  is  substantially  less  diffusive,  only  slightly  more  expen¬ 
sive,  and  quite  robust.  The  Roe  scheme  is  the  least  diffusive,  costs  about 
50%  more  than  the  Boltzmann  scheme,  but  is  not  as  robust  as  the  HLLL 
scheme.  Figures  14  and  15  show  detailed  results  for  the  Roe  scheme. 

In  problems  in  which  the  magnetic  field  is  particularly  strong,  the  total 
energy  can  be  dominated  by  the  B  •  B  term.  Often  in  these  cases,  however, 
the  strong  magnetic  field  is  very  close  to  a  known  solution,  for  example,  a 
dipole  field.  Writing  the  governing  in  terms  of  the  perturbation  from  a  known 
field  Bo  eliminates  the  Bq  •  Bo  terms  in  the  equations,  thereby  eliminating 
some  numerical  difficulties.  It  is  not  unusual  for  Bo  to  be  several  orders 
of  magnitude  larger  than  the  perturbation  field;  this  subtraction  technique 
substantially  increases  the  accuracy  with  which  the  perturbation  field  is  re¬ 
solved.  It  is  important,  however,  to  do  this  subtraction  intelligently.  The 
resulting  Riemann  solver,  in  which  the  unknowns  are  based  on  the  perturba¬ 
tion  magnetic  field,  has  been  developed,  tested,  and  written  up  for  submission 
to  a  journal  (Reference  1).  The  treatment  that  allows  this  approach  to  be 
parallelized,  achieving  linear  speed-ups  even  on  hundreds  of  processors,  was 
reported  at  the  Supercomputing  conference  (Reference  2). 

The  basic  geometry  generation  capabilities  of  the  code  developed  by 
Charlton  (Reference  3)  have  been  extended  to  allow  parametric  definition 
of  arc-constrictor  geometries.  This  tool  allows  generation  of  these  geometries 
with  almost  any  combination  of  dimensions  in  a  matter  of  minutes.  Initial 
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grids  are  generated  in  a  matter  of  hours  using  automated,  adaptive  Cartesian 
methods,  and  solution-based  adaption  allows  the  code  to  increase  resolution 
around  flow  regions  of  interest  such  as  the  arc  during  flow  solution.  This 
allows  application  of  the  code  to  a  large  range  of  problems  involving  MHD 
and  complex  geometries;  arcjets,  magnetoplasmadynamic(MPD)  thrusters, 
and  fusion  propulsion  systems  are  just  a  few  examples. 
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Constricto 


Figure  3:  Triple  Probe  Circuit  Used  for  Near-Field  Plume  Measruments. 
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Figure  4:  High-Frequency  Voltage  and  Current  Oscilloscope  Traces. 
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Figure  5:  High-Frequency  Voltage  and  Current  Oscilloscope  Traces* 
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Figure  7:  High-Frequency  Voltage  and  Current  Oscilloscope  Traces. 
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Figure  8:  Electron  Number  Density  Axial  Profile. 
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Figure  9:  Electron  Number  Density  and  Temperature  Profile. 
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